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The dynamics of TV point vortices in a fluid is described by the Helmholtz-Kirchhoff (HK) equa- 
tions which lead to a completely integrable Hamiltonian system for TV = 2 or 3 but chaotic dynamics 
for TV > 3. Here we consider a generalization of the HK equations to describe the dynamics of mag- 
netic vortices within a collective-coordinate approximation. In particular, we analyze in detail the 
dynamics of a system of three magnetic vortices by a suitable generalization of the solution for three 
point vortices in an ordinary fluid obtained by Grobli more than a century ago. The significance of 
our results for the dynamics of ferromagnetic elements is briefly discussed. 



I. INTRODUCTION 

Ferromagnetic media are described by the density of magnetic moment or magnetization 
m = m(r,t) which is a vector field of constant length that satisfies the highly nonlinear Landau- 
Lifshitz (LL) equation. Magnetic bubbles and vortices are special solutions of the LL equation 
characterized by a topological invariant sometimes called the Pontryagin index or skyrmion num- 
ber. Upon suitable normalization, the skyrmion number is integer for magnetic bubbles occuring 
in magnetic films with perpendicular easy-axis anisotropy [l[ and half integer for magnetic vor- 
tices that may be realized in the case of magnetic films with easy-plane anistotropy [2J • 

Most of the early studies were carried out on magnetic films of infinite extent where vortices 
and bubbles occur as excited states above a uniform ground state. More recently it has been 
realized that the ground state of a disc-shaped magnetic element, with a diameter of a few 
hundred nanometers, is itself a vortex configuration. Hence the vortex is a nontrivial magnetic 
state that can be spontaneously created on magnetic elements of finite extent. Such a possibility 
is of obvious significance for device applications and explains the current interest on this subject 

Hi. 

The dynamics of magnetic vortices is greatly affected by the underlying topological structure [fj, 
7] . This is described by a topological vorticity which gives the skyrmion number upon integration 
over the plane of the magnetic film. As a result, the essential features of magnetic vortex 
dynamics are similar to those displayed by ordinary vortices in fluid dynamics. A single vortex 
or antivortex is spontaneously pinned and thus cannot move freely within the ferromagnetic 
medium. However, vortex motion relative to the medium is possible in the presence of other 
vortices, or externally applied magnetic field gradients, and displays characteristics analogous 
to those of two-dimensional (2D) motion of electric charges in a uniform magnetic field. In 
particular, two like vortices orbit around each other while a vortex-antivortex pair undergoes 
Kelvin motion. A recent review of two-vortex magnetic systems together with some progress in 
the study of three- vortex dynamics may be found in Refs. [1, Q. 

In the present paper we study the dynamics of a system consisting of TV interacting magnetic 
vortices within a collective-coordinate approximation which evades the complexities of the com- 
plete LL equation but retains important dynamical features due to the underlying topology. The 
relevant equations are a generalization of the HK equations encountered in the study of TV point 
vortices in ordinary fluids [lol [Tlj . The latter equations were solved explicitly by Helmholtz 
himself for the two-vortex (TV = 2) system, while Kirchhoff concluded that the three-vortex 
(TV = 3) system is completely integrable. But a complete solution for TV = 3 was given in the 
1877 dissertation of Grobli [n| recently reviewed by Aref et al [l3[ (see also For TV > 3, 

the system is not completely integrable and thus leads to chaotic dynamics. 
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In Section II, we state the appropriate modifications of the HK equations to account for N 
interacting magnetic vortices, derive the corresponding conservation laws, and provide a hrst 
demonstration with an explicit solution for two-vortex systems. The generalization of Grobli's 
solution to three magnetic vortices is given in Section III and is illustrated in detail in Section 
IV for important special cases such as scattering of a vortex-antivortex pair off a target vortex 
initially at rest, three-vortex collapse, bounded three-vortex motion, etc. Our main conclusions 
are summarized in Section V. 



II. GENERALIZATION OF THE HELMHOLTZ-KIRCHHOFF EQUATIONS 

We consider a system of N magnetic vortices labeled by a = 1,2, ...N. Each vortex is 
characterized by a pair of indices (k q , A q ) where K a is the vortex number and A Q the polarity, 
which take the values n a = ±1 and A a = ±1 in any combination. We also define the skyrmion 
number 

S a = K a X a (1) 

which also takes the distinct values s a = ±1 (in any combination for varying a = 1,2,.. . N) 
and differs from the standard definition of the Pontryagin index by an overall normalization 

{N a = -± Sa ) !!. 

Now, in a collective-coordinate approximation [15|, [l6|, [l7| , the position of each vortex is labeled 
by a vector r a = (x a , y a ) in the 2D plane and the energy functional is given by 

E = - K a K/3 In \r a - rp\ (2) 

Q</3 

in suitable units. Note that the energy is a function of the vortex number n a but not the polarities 
A Q . The equations of motion are then written in Hamiltonian form as 



dx a dE dy a dE 

a dt dy a ' a dt dx a 

or, more explicitly, as 



1,2,. ..N, (3) 



dx a \ - Va-y/3 , dy a \ - x a - Xp 

Aa "~£ K/3 K^F' |>F^F- (4) 

These differ from the HK equations for N point vortices in an ordinary fluid by the presence of 
the polarities Aq, (or the skyrmion numbers s a ) in addition to the vortex numbers K a . Equations 
(|4]) reduce to the HK equations in the special limit where all polarities are taken to be equal 
(Ai = A2 = • • • = Ajv). Also note that the vortex number n a may assume any continuous value 
for a fluid vortex but takes the quantized values ±1 for a magnetic vortex. 

In addition to the conserved energy of Eq. @, the Hamiltonian system ((3J) or Q possesses 
conservation laws which originate in translational and rotational invariance; i.e., the analog of 
linear momentum P = (P x , P y ) with 

Px ^ $a y<y-> Py ^ Sa ^ott (5) 

a a 

and angular momentum 

L = \Y,^^l + yl) (6) 
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both defined in convenient units. Actually, the physical significance of these conservation laws 
depends crucially on whether or not the total skyrmion number 

S = ^ Sg = ^ Kg A a (7) 

a a 

vanishes. The linear momentum ([5]) assumes its customary meaning only when S = 0. For 
S 0, one may instead define the conserved guiding center R = [R x , R y ) with 

Re — ~q ^ 5 q^qi Ry ^ ^ ^aj/q ; (8) 

a a 

which is a measure of position rather than linear momentum. The inherent transmutation of 
momentum into position occurs also within the complete LL equation and has been the subject of 
several investigations [f| 0, • To complete the discussion of conservation laws we note that the 
linear momentum ([5]) and the angular momentum (JSJ) may be combined to yield the derivative 
conserved quantity 

\ J2s aSfi {r a -Tpf (9) 

a 13 

which depends only on the relative distances \r a — rp\ and will play an important role in the 
following. 

This section is concluded with a simple application of the general formalism to the case of 
two- vortex systems (N — 2). Equations ^ then read 

dxi yi-V2 , dx 2 _ V2-yi 

dyi xi-x 2 , dy 2 x%-x x 
Al ^ = K2 ^^' A2 " Kl ^^' (10) 

where d 2 = (x 2 — xi) 2 + (y 2 — y\) 2 is the squared relative distance between the two vortices which 
is conserved by virtue of the conservation of the energy E of Eq. @ restricted to N = 2 for any 
choice of vortex numbers and polarities. But the details of the motion depend on the specific 
choice of k's and A's as demonstrated below with two representative examples. 
First, we consider a vortex-antivortex pair with equal polarities: 

(M,A 1 ) = (1,1), (/sa.Aa) = (-1,1), (11) 

hence opposite skyrmion numbers s\ — 1, s 2 = — 1 and total skyrmion number S = s x + s 2 = 0. 
A straightforward solution of Eqs. (fT0)l then leads to Kelvin motion where the vortex and the 
antivortex move in formation at constant relative distance d, along parallel trajectories that are 
perpendicular to the line connecting the vortex and the antivortex, with speed 

V = \ (12) 

and linear momentum ([5]) that is parallel to the velocity. 

In contrast, a vortex-antivortex pair with opposite polarities: 

( Kl ,A 1 ) = (-l,l), (k 2 ,A 2 ) = (1,-1), (13) 

hence s\ = s 2 = —1 and nonvanishing total skyrmion number S = s\ + s 2 = —2, undergoes 
rotational motion at constant diameter d, around a fixed guiding center ([5]) which may be taken 
to be the origin of the coordinate system {x\ + x 2 = = yi + y 2 ), with angular frequency 

-J, (14) 
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A detailed study of Kelvin and rotational motion was carried out within the complete LL equation 
in Refs. [l8[ and [T^|, respectively, the results of which were recently reviewed in Ref. [8|. In 
particular, their results confirm the validity of Eqs. (fT2)l and (fl4|) for sufficiently large relative 
distance d, as expected for the collective-coordinate approximation adopted in the present paper. 

Finally we note the the motion of a two-vortex system is Kelvin-like when the total skyrmion 
number vanishes (S = s\ + s 2 = KiX\ + K2X2 = 0) and rotational otherwise (S = si + s 2 = 
«iAi + K2X2 = ±2) irrespectively of the detailed configuration of the vortex numbers K\, K2 and 
polarities Ai, A2. 



III. THREE MAGNETIC VORTICES 



Consider a system of three interacting magnetic vortices with vortex numbers and polarities 
given by the three pairs of indices (ki, Ai), (k 2 , A2), (k 3) A3). The three vortices form a triangle 
whose vertices are located at r\ — (£1,3/1), r 2 = (£2,2/2), ^3 = {x 3 ,y 3 ) and move about the 2D 
plane according to Eqs. written here explicitly as 

dxi yi - y 2 2/3-2/1 x dyi _ X\ - x 2 x 3 - x\ 

dx 2 2/2 - 2/3 , 2/1-2/2 . dy 2 x 2 - x 3 Xi - x 2 . . 

A 2 = - K 3 C 2 + «1 ^2 » 2 ~dT^ K3 ~Cf Kl g 2 ■ ( 15 ) 

, dx 3 2/3 - 2/1 , 2/2-2/3 . dy 3 _ x 3 -x x x 2 - x 3 
A 3 — — - ~ K i — 7^2 r K z — Ta — ' A 3 -jr - hi — —2 n 2 - 



where 



dt x ci C( 1 " dt 1 CI ' C( 



Ci = |r2-r 3 |, C 2 = |r 3 -ri|, C 3 = \r 1 -r 2 \ (16) 



are the lengths of the triangle sides. 

Our task is to solve the initial- value problem defined by Eqs. (j 1 5[) with initial conditions 
furnished by the t = configuration of the three vortices. Needless to say, a numerical solution 
is straightforward and has indeed been used to confirm our main analytical results obtained in 
the following by a generalization of Grobli's original solution for ordinary fluid vortices [l2l ]. 

The key observation is that the scalar quantities C\ , C2 and C 3 satisfy a closed system of 
equations of their own, which may be employed to determine the time evolution of the shape 
and size of the vortex triangle, irrespectively of its relative orientation and motion in the plane. 
Specifically, as a consequence of Eqs. ([T5)| . 

d {CD = A Kl A 



1 


1 


A 3 C| 


A 2 Cf 


1 


1 


Ai Cf 


A3C? 


1 


1 


A2C? 


AiC 2 



dt 

|(C 2 2 )=4^(— ) (17) 
j t (C!)=4 K3 A 

Here A = i/\A\ where v = 1 or —1 when the vortices (1,2,3) appear in the counterclockwise or 
clockwise sense in the plane, and \A\ is the area of the vortex triangle which may be expressed 
entirely in terms of C\ , C2 and C 3 : 

\A\ = [C(C-C 1 )(C-C 2 )(C-C 3 )] 1 I 2 , C=^(C 1 + C 2 + C 3 ). (18) 
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Eqs. (|17[) constistute a reduced dynamical system that is completely integrable because it pos- 
sesses two conservation laws; namely 

B = c 1/Kl c 1/K2 c 1/K3 r = S2S3C i + s * SlC l + sis 2 Cl 

1 2 3 ' Si + S 2 +S 3 

where the conservation of B follows from the conservation of the energy of Eq. and V from 
the conservation of the quantity defined in Eq. ([9]) , both restricted to the three- vortex system 
and normalized in a manner suitable for subsequent analysis. The conservation of B and V may 
be confirmed directly from Eqs. (TIT)) . Thus, having determined B and T from the initial (t = 0) 
values of Ci, C 2 and C 3 , one may use Eqs. (|19[l to eliminate two of these variables, say, C 2 and 
C3 in favor of C\ and thereby reduce Eqs. (|17[) to a single equation for C\ that may be integrated 
by a simple quadrature. 

In order to complete the solution we must also determine the relative orientation and motion 
of the vortex triangle in the plane and thus solve for the actual vortex trajectories. The latter 
may be conveniently described in polar coordinates (r a , 9 a ), a = 1, 2,3. Using these coordinates 
Hamilton's equations © read 

dr^_dE_ d9« _dE_ 

SaTa dt ~ de a SaTa dt _ dr a [ U) 

We write explicitly the three equations for the polar angles: 

dOi n - r 2 cos(0i - $2) 

Ai n — = n 2 h K3 

dt 

d9 2 r 2 - r 3 cos(6> 2 - 63) , • .,. 1 

^r 2 — = K 3 + «i — 2 . <- >Jl 

d9 3 r 3 - n cos(0 3 - 9 X ) 
dt C2 

A crucial step for simplifying the latter equations is to set the guiding center (JSJ) at the origin, 
i.e., R x = = R y . In polar coordinates this constraint reads: 

sir 1 cos#i + s 2 r 2 cos 6 2 + s 3 r 3 cos #3 = 

Sxfi sin 0i + S2T2 sin #2 + s 3 r 3 sin 9 3 = 0. (22) 
Simple combinations of Eqs. (f2"2"]) yield 

2s 2 s 3 r 2 r 3 cos(6> 2 - 3 ) = s\r\ - s\r\ - s\r\ 

2s 3 S! r 3 ri cos(6» 3 - 6»i) = -sjrj + s\r\ - s\r\ (23) 
2sis 2 rir 2 cos(#i - 2 ) = -s\r\ - s\r\ + Sgfg. 
Further, we write the sides of the triangle using polar coordinates: 

C\ = r\ + rt - 2r 2 r 3 cos(0 2 - 6 3 ) 

Cf = r l + r l - Irzri cos(6» 3 - 9 X ) (24) 

C 3 = r\ + r\ - 2nr 2 cos(#i - 2 ), 

and we can now use Eqs. (|2"3"|) and (fM]) to obtain explicit relations between the radial coordinates 
r a and the sides of the triangle C a : 

s 2 s 3 Cl = (s 2 + s 3 )T - si(si + s 2 + s 3 )rj 

S3S1 G\ = (s 3 + si) r - s 2 (si + s 2 + s 3 ) r£ (25) 
sis 2 Cf = (si + s 2 ) T - s 3 (si + s 2 + S3) fg, 



ri 


- r 3 cos(6» 3 - 


61) 




c\ 




r 2 


— n cos(#i — 






CI 




''3 


— r 2 cos(# 2 — 


e 3 ) 


c? 
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where T is the conserved quantity already defined in Eq. (|19p . 

We are now ready to write Eqs. ([2"Tj) in a form where the radial coordinates r a are eliminated 
and only the variables C a and 6 a appear. We use Eqs. ((24)) to eliminate the cosines and Eqs. ([25]) 
to express the radial coordinates in terms of the C a 's in Eqs. (|2"Tj) . to obtain 



where it is understood that the r Q 's are given through Eqs. (f2"S")) in terms of the C Q 's. Explicit 
applications of the preceding general results are given in Section IV. 



For an explicit demonstration we consider a three-vortex system with vortex numbers and 
polarities given by 



and corresponding skyrmion numbers s\ = 1, s 2 = —1 = S3. The pair (1,2) is the vortex 
-antivortex pair of Eq. (jTTJ) which would undergo Kelvin motion in the absence of vortex 3. In 
contrast, the pair (2, 3) would undergo rotational motion in the absence of vortex 1, as discussed 
in Sec. II following Eq. (Ti"3"|) . Thus the motion of the three vortices is expected to be a combination 
of Kelvin and rotational motion. Since the total skyrmion number S = Si + s 2 + S3 = — 1 is 
different from zero, the guiding center for the three- vortex system defined from Eq. l[8"]). i.e., 



2A 1 A 2 A 3 (si + S2 + sz)r\ClCl -± = X 3 C 2 [s 2 s 3 (-C 2 + C* 2 + C 2 ) + 2s\C 2 \ 

+ A 2 C 2 [s 2 s 3 (-Cf + C\ + Cl) + 2 S 2 C|] 

2A 1 A 2 A 3 ( Sl +s 2 + s^r 2 2 ClCl ^ = X^l [« 3 *i(C? - C\ + C 3 2 ) + 2s 2 3 Cf] 

+ A 3 Cf [ S3S i(C 2 - Cl + Cl) + 2s\Cj] 

2A!A 2 A 3 ( Sl +s 2 + s 3 )rjCjC 2 ^ = A 2 C 2 [ Sl s 2 (C 2 + C 2 - C 2 ) + 2s 2 C 2 ] 

+ X1C 2 [s lS2 (C 2 + C 2 - C 2 ) + 2s 2 C 2 }, 



(26) 



IV. AN INTERESTING THREE- VORTEX SYSTEM 



(ki,Ai) = (1,1), (k 2 ,A 2 ) = (-1,1) > (k 3 , A3) = (1,-1) 



(27) 



R x = -xi + x 2 + x 3 , 



Ry 



Vi +2/2 + 2/3, 



(28) 



is conserved and thus remains fixed during the motion. 

We now return to the reduced system (jTTJ) applied for the specific choice (12~T1) : 




(29) 



while the conservation laws defined from Eq. (|19[) are given by 



B = 



Ci C 3 

c 2 



r = -cl + cl + c 3 2 . 



(30) 
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One may then express C2 and C3 in terms of C\ as 

2 _ C\ + r 2 ^,2 _ Cl + T R 2 /01 \ 

and the area of the triangle defined earlier in Eq. (fT5|) as 



4A = v^ClCl - F» = V 4B2C i 2 ( C i 2 + r ) 2 - r2 ( C i 2 + S2 ) 2 ' ( 32 ) 

Then the system of Eqs. (|2"5)) reduces to the single equation 

|(c?) = -^pr (cf + \/ 4i?2c i( c i +r ) 2 - r2 ( c i +B2 ) 2 ' ( 33 ) 

which may readily be integrated to yield C\ — C\ (t) , as well as C2 — C2 (t) and C3 = C3 (i) from 
Eq. (131]) . for any specific values of the conserved quantities B and T calculated from the initial 
three-vortex configuration. 

In order to complete the solution we must also determine the relative orientation and motion 
of the vortex triangle in the plane and thus solve for the actual vortex trajectories. We will 
follow the derivation given in Section III and quote only important results that will allow us to 
provide a detailed illustration of the three- vortex motion for the special case of vortex numbers 
and polarities given by Eq. (|2"T)) . 

It is convenient to specify a reference frame such that its origin coincides with the conserved 
guiding center of Eq. (|28| : 

-xx + x 2 + x 3 = 0, -yi +2/2+2/3 = 0. (34) 

In this frame the radial distances r\ , r-i and r$ are related to the triangle sides C\ , C2 and C3 by 
Eq. (p?5l) applied for Si = 1, S2 = — 1 = S3: 

t\ =2(Cl + Cl)-Cl, rl = Cl, r 2 3 =Cl (35) 

while the angular variables satisfy Eqs. (|26p now reduced to 

d8 1 r 



dt 2{C\ + 2r) 

^1 = —( — 

dt ~ 2C| \Cl 

d6 3 r / 1 









C3 2 ) 




1 




+ cl 




1 


c?)" 





(36) 



To summarize, for the specific choice of vortex numbers and polarities given by Eq. (|27j) . the 
radial variables n , ri and r3 are expressed in terms of C\ , C2 and C3 from Eqs. (|35p and ultimately 
in terms of Ci alone using the conservation laws as in Eqs. (|3ip ; while the time evolution of C\ 
is determined from Eq. ([3"3"]l by a straightforward integration. Similarly, we may use Eqs. ([3~lj) 
to express the right-hand sides of Eq. |36|) in terms of C\ and thus the time evolution of the 
angular variables 9\, 62 and #3 is also reduced to simple quadratures. 

Without loss of generality, the initial three-vortex configuration is shown in Fig. [TJ In par- 
ticular, b is the initial length of the base (1,2) of the triangle and h the corresponding height. 
The impact parameter a is defined as the distance of vortex 2 from the y-axis and is taken to be 
negative when vortex 2 is located to the left of the y-axis and positive otherwise. Also note that 
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FIG. 1: Initial configuration of three magnetic vortices (1,2,3) with vortex numbers and polarities given 
by Eq. (I27|l . The origin of the coordinate frame is taken to coincide with the conserved guiding center 
of the three-vortex system. The impact parameter a is taken to be negative when vortex 2 lies to the 
left of the y-axis and positive otherwise. 



vortex 3 is initially located on the negative x-axis at a distance 6 from the origin, thus enforcing 
Eqs. (|3"4"1) . Hence, the initial values of the triangle sides are C\ — h 2 + (a + 6) 2 , C 2 = h 2 + a 2 , 
and C 2 = b 2 , and the conserved quantities of Eq. ([50]) are given by 

B>= h2 + { " + 2 b)2 #, T = -2ab, (37) 
h z + a 1 

where T is positive for left impact (a < 0) and negative for right impact (a > 0). This completes 
the description of the general procedure which is explicitly implemented in the continuation of 
this section using examples of increasing complexity. 



A. Special case 

First we employ the general formalism developed in the present paper to recover a special 
solution already given in the Appendix of our earlier paper [9(. The initial configuration is 
taken to be an isosceles triangle defined by choosing the impact parameter a = —6/2 in Fig. [1] 
The vortex-antivortex pair (1, 2) would then move upward with initial velocity V ~ 1/6 (Kelvin 
motion) and eventually collide against the target vortex 3. After collision the vortex-antivortex 
pair moves off to infinity at some scattering angle while the target vortex comes to rest at a new 
location. Our aim in the following is to calculate the scattering angle as well as the final position 
of the target vortex. 

Now, inserting a = —6/2 in Eq. ([57)) . we find that the values of the conserved quantities are 
given by 

B 2 = b 2 =T (38) 

and are independent of the initial triangle height h. Then Eqs. (f3"Tj) read 

C l (t)=C 2 {t), C 3 (t)=b (39) 

which establish that the triangle remains isosceles (Ci = C2) and the length of its base remains 
constant (C 3 = 6) at all times t > 0. But C\ itself undergoes a nontrivial time evolution governed 
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by Eq. HMD: 



dt^ u \C\ ' b 2 

A more convenient parametrization is obtained by 



d (C!) = -vb(± + ±)jAC>-l. ,40, 



C 1 = WiP + I&2 = C2) C 3 = 6, (41) 

where H = Hit) is the instantaneous height of the triangle. Eq. (|4Hll then reads 

dH_ V R- + W 

dt b H 2 + \b 2 V ' 

and must be solved with initial condition H(t = 0) = h. 

During the initial stages of the evolution the three vortices (1,2,3) appear in the plane in a 
counterclockwise sense and thus v = 1. Eq. (|42p is then integrated to yield 

H 4 fH\ t -t 

___ arctan /_j = _ 



arctan) — ) = < t < t (43) 



h 4 / h 

C6-~5 arCtan U 



As t approaches to, height H approaches zero and the triangle reduces to a straight line element 
(1,3,2) at t = to with vortex 3 located at the center of the line element (1,2). Also recalling 
that the length of the line element (1,2) remains constant (equal to b) we conclude that to is 
the instance of the closest encounter of the three vortices. For t > to, vortex 3 emerges from the 
other side and thus Eq. (|4"2")) must be integrated with v = — 1 and initial condition H(t = to) = 0: 

H 4 fH\ t-t 

ti arctan — = — — 



— - - arctan ( — ) = —r^, t > t . (44) 



In the far future (t — > oo) H — > t/& which is consistent with a vortex-antivortex pair (1,2) moving 
away from the target vortex 3 with asymptotic speed V = l/b, in agreement with Eq. Q12p 

The preceding results already suggest a process in which the vortex-antivortex pair (1,2) 
initially approaches the target vortex 3 but eventually moves off to infinity at a scattering angle 
that remains to be calculated. One must also ascertain the fate of the target vortex 3 well after 
collision. Both of these questions are answered below by applying Eqs. (|35[) and (|36p to the 
special case treated in this subsection. 

First, we note that Eqs. (|35[) suggest that the three vortices (1,2,3) together with the origin of 
the coordinate system O (i.e., the guiding center of the three- vortex system) form a parallelogram 
1203 at all times during the collision. In particular, the third equation establishes that 
7'3 = C3 = b is time independent, or, vortex 3 moves on a circle whose center coincides with 
the guiding center and its radius is a constant b equal to the constant length of the base of the 
triangle. Furthermore, the polar angle #3 changes during collision by the same amount as the 
base of the triangle (1,2) rotates in the 2D plane. Therefore, the total change A83 is equal to 
the scattering angle of the vortex-antivortex pair (1,2) and also specifies the final location of the 
target vortex 3. 

In order to actually calculate A6*3 in the present special case, we apply Eq. ([55]) with C\ = 

d63_J__ 1 

dt c\ m + \b^ { ' 
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FIG. 2: Time evolution of C\,C% and C3 for an initial three- vortex triangle given by Fig. [T]with a — 
0, b = 1, h = 5. Note that all three lengths simultaneously collapse to zero at a finite time to calculated 
from Eq. (53]). 



where we may insert H = H{t) from Eqs. (|4"5|) and (|4"4"]) to obtain 6*3 = 6^{t) by a further 
integration. In fact, a more direct calculation of the total scattering angle A83 is obtained by 
combining Eqs. (JUJ) and (|4l)]) to write 



d6>i 6 1 



and 



dH v H 2 + |6 2 ' 



(46) 



M3 = Vo lPTW + h k WTW> (47) 

where the two terms may be interpreted as the scattering angles before and after collision and 
correspond to the two branches of the solution H = H(t) given by Eqs. (|43[) and (J44]) . In 
particular, for pure scattering where the vortex-antivortex pair originates very far from the 
target vortex {h — * 00): 

A6 3 = 2 b r 2 dH 5 =^, (48) 

which reproduces the special result given in the Appendix of Ref. 0] . A detailed illustration of 
the actual vortex trajectories may be found in Fig. Al of the above reference (using a different 
convention for the origin of the coordinate system) and in Fig. [5] of the present paper (where the 
origin of the coordinate system coincides with the guiding center). 
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B. Three- vortex collapse 



At first sight, the special example treated in the preceding subsection would seem to cor- 
respond to a head-on collision of a vortex-antivortex pair off a target vortex. However, such 
an interpretation seems to be unjustified in view of the peculiar scattering angle calculated in 
Eq. (|48p . Indeed, an example that is closer to head-on collision is obtained by choosing vanishing 
impact parameter (a = 0) in the initial vortex configuration shown in Fig. [TJ 

The initial vortex triangle is then orthogonal, with initial values of the length of its sides 
Ci(t = 0) = \/h 2 + b 2 , C 2 (t = 0) = h and C 3 {t = 0) = b. The two conserved quantities of 
Eq. (|3"0)) now read 



5 = ^ = 70?+^, r = -C\ + C\ + C\ = 0. (49) 

A notable consequence of the conservation of T = is that the vortex triangle remains orthogonal 
at all times, even though the length of its sides Ci, C2, and C3, as well as its relative orientation 
in the 2D plane, undergo a nontrivial time evolution. 

In particular, the evolution of C\ is governed by Eq. (|3"3")l applied with T = and B given by 
Eq. (011): 

dC 1 Cl + B 2 

~df = ~^cT' m 

while C2, and C 3 obtained from Eq. (|3ip are now given by 



C2 -7W^' ^-^TW- (51) 



Eq. (|5D|) may easily be integrated to yield 



§- arCtan (§)= V' (52) 



where the integration constant £o is calculated from the initial condition C\ (t = 0) = Vh 2 + b 2 , 
or 



t = ^ 2 (h 2 + b 2 ) 



h (h 

arctan — 

b \b 



(53) 



The physical significance of £o becomes apparent by applying Eq. (|52p in the limit t — > to: 

d(t^to)~ [3B(t Q -t)}^ 3 , (54) 
and Eqs. (|51[) in the same limit to yield 

C 2 (t^t )~ ±-{3B(t -t)} 2 /\ C 3 (t^t ) ~ [3B(t -i)] 1/3 . (55) 

Therefore, to calculated explicitly from Eq. |53|) is the instance at which all three C\ , Ci and 
C3 vanish simultaneously, or, the vortex triangle collapses to a point. Also taking into account 
Eq. ([55)) , we conclude that the vortex triangle collapses to the guiding center (ri = r 2 = r 3 = 0) 
of the three- vortex system in a finite time interval to. The exact time dependence of C\, C 2 and 
C3 calculated from Eqs. ([52")) and (|5T|) over the entire time interval < t < to is depicted in 
Fig.H 
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FIG. 3: Trajectories of the three vortices 1,2 and 3 emanating form an initial three-vortex configuration 
given by Fig. [T] with a = 0, b = 1, h = 5. Note that the three trajectories merge at the origin of the 
coordinate frame which is taken to coincide with the guiding center of the three- vortex system. 



Complete information about the vortex trajectories may be obtained by further utilizing the 
results of Eqs. ([55]) and or by a direct numerical integration of Eqs. (jTSJ) . The actual 
trajectories are illustrated in Fig. [3] which explicitly demonstrates the predicted three-vortex 
collapse. A notable feature of Fig. [3] is that vortex 1 moves on a straight line that connects its 
original position with the guiding center of the three-vortex system. 

Although this special case of vanishing impact parameter (a = 0) is clearly singular from the 
point of view of scattering theory, it is still possible to calculate the scattering angle A#3 before 
collapse. First, we apply the third of Eqs. Q36p with T = 0: 

~A=Cl (56) 



and combine this result with Eq. (f5T}| to obtain 



Bdd 
B 2 + C 2 



dO* = -vTT-k- (57) 



The scattering angle before collapse is then obtained by integrating Eq. ([57]) between the values 
Ct(t = 0) = Vh 2 + b 2 and d(t = 0) = 0: 



Bdd (h 



A0 3 |before = J ^ 2 ~ ^ = arctan \ t J ■ (58) 

In the limit where the vortex-antivortex pair originates at a large distance h — > oo from the 
target vortex, we find A#3|bcfore = and, by convention, A03| total = 2(7r/2) = 7r. This result 
agrees with the interpretation of this special case as a head-on collision whose singular nature 
will also become apparent in the discussion of the following subsection. 
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C. Scattering in the general case 



Our aim in this subsection is to determine the scattering angle as a function of impact pa- 
rameter a. The general strategy is suggested by the special examples treated in the preceding 
two subsections. Thus we return to the initial configuration depicted in Fig. [TJ now applied for 
arbitraty a, which yields values for the conserved quantities B and T as given in Eq. (|37p . In 
particular, we shall be interested in the limit h — > oo: 

B 2 = b 2 , r = -2ab, (59) 

which corresponds to pure (asymptotic) scattering where the vortex-antivortex pair is initially 
located very far from the target vortex and also diverges to infinity after collision. 

In order to calculate the scattering angle we return to Eq. (f3T>)) and express its right-hand side 
entirely in terms of C\ using Eq. (|3"Tj) to eliminate C2 and C3. We then combine the resulting 
equation with Eq. (|33|) to write 



d0 3 __ 1 2cf+r + c'(+B* a 



dCf v ^/AB 2 C 2 (C 2 + T) 2 - T 2 (C 2 + B 2 ) 2 ' 
Simple inspection of this equation suggests using scaled variables u and 7 from 



(60) 



C? r 2a 

7 = ^ = ~t (61) 



to obtain 



, n 7(1-7) 1 M+7 

ac/3 _ 1 2(u+7) ' u+1 



du v ^/Au(u + j) 2 - 7 2 (w + l) 2 



(62) 



Based on this equation the total scattering angle is calculated as a function of the parameter 7 
alone by a procedure similar to that employed in the special example treated in subsection A: 

f oc 7(1-7) , u+7 

A9 3 = 2 2(u+l) u+1 du (63) 

Juo v/4 M ( u + 7) 2 -7 2 ( u + l) 2 

where the factor of two and the upper limit of integration are dictated by the fact that the vortex- 
antivortex pair is located very far from the target vortex well before and after the collision. The 
lower limit uq is a specific value of u = C 2 /B 2 for which the area of the vortex triangle vanishes. 
Therefore, uq must be chosen among the three roots of the cubic equation 

4?i(?i + 7 ) 2 -7 2 (u + l) 2 = 0. (64) 

It is thus important to examine the behavior of the roots as functions of the parameter 7. 

For 7 < there exists a real and positive root u\ and two complex roots ui and 113 such that 
U2 = W3. In this region the lower limit in the integral of Eq. (|63[) must be chosen as uq = u\. 
Now, 7 = is an exceptional point in that all three roots vanish (u\ = U2 = u% = 0) and thus 
the lower limit must be chosen as up = 0. Applying Eq. ([6"3"]) for 7 = and uq — yields 

A9 3 = Z r =« (65) 

Jo (w+l)V u 

in agreement with the conclusion of Subsection B. Indeed, 7 = corresponds to vanishing impact 
parameter (a = 0) which is the special case discussed in Subsection B (three- vortex collapse). 
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Next we consider the region < 7 < 1 where the cubic equation again possesses a real and 



positive root u\ and two complex roots U2 



Thus the lower limit must again be chosen 



as uq = u\. Now, 7 = 1 is another exceptional point in that all three roots become real 
[u\ = j,U2 = —1 = U3) and the lower limit must be chosen to coincide with the positive root 
uq = u\ = 4. Applying Eq. (j6"5)) for 7 = 1 and uq — \ yields 



A0 3 = 2 



du 



1 



2tt 

71' 



(66) 



which coincides with the result of Eq. (14"5|) . This is not surprising because 7 = —2a/b= 1 leads 
to an impact parameter a = —6/2 which is indeed the special case discussed in Subsection A. 

For 7 > 1, the cubic equation again possesses a real and positive root u\ and two complex 
roots U2 = W3, thus we must choose uq = u\. This generic picture continues to hold until 7 
reaches the critical value 



7c 



27 

y 



(67) 



where all three roots become real {u\ — 0.5625, U2 = 9 = M3) and remain real and positive for 
7 > 7 C . In this range of 7, the three roots may be ordered as u\ < 112 < U3 and the lower limit 
of integration in Eq. (|63[) must now be chosen as uq — U3. 
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FIG. 4: Scattering angle during collision of a vortex-antivortex (Kelvin) pair (1,2) with a target vortex 
3, as a function of impact parameter a. Note the crossover behavior at the critical values a/b = and 
a/b — —27/4 discussed in the text. The results depicted in this figure should be read modulo 2ir. 



In all cases the denominator in Eq. (|fc>3[) does not possess singularities in the domain of inte- 
gration [wo,oo) except for an integrable square root singularity at the lower end. It can also be 
shown that the numerator in Eq. (|6"5)) does not exhibit singularities in the integration domain. 
Thus we have provided a complete prescription for the calculation of the scattering angle, which 
entails locating the appropriate root uq = uq (7) of the cubic equation and an elementary numer- 
ical integration of the integral in Eq. (|6"5)) for any value of 7 = —2a/b. The calculated scattering 
angle as a function of impact parameter a measured in units of b is depicted in Fig. [4] 

Thus the scattering angle displays critical (crossover) behavior at two characteristic values of 
the impact parameter; namely, a/b = 0, which corresponds to the three- vortex collapse discussed 
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FIG. 5: Vortex trajectories for six characteristic values of the impact parameter a. In all cases a vortex- 
antivortex (Kelvin) pair (1,2) collides with a target vortex 3 and is deflected by a scattering angle in 
agreement with the results of Fig. [4] The target vortex 3 is at rest well before collision and also comes 
to rest well after collision, as expected for an isolated vortex which is always spontaneously pinned. 



In Fig. [5] we display the actual vortex trajectories calculated numerically from Eqs. (|15[) for a 
set of six characteristic values of a/b. The change of pattern becomes clear as one crosses the 
two critical values a/b — and a/b — —27/4. The nature of the critical parameter 7 C of Eq. ([57]) 
is further illuminated in the following subsection. 
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D. Bounded three-vortex motion 



The analysis of the cubic equation (|64[) revealed the importance of the critical parameter 
7 C = 27/2 of Eq. ([67)) above which all three roots are real (and positive). We may then order 
the roots as u\ < u 2 < u 3 and note that the argument under the square root in, say, Eq. (|62l) is 
positive when either u > u 3 or u\ < u < u 2 . The former case, 

27 

7 > 7c = y, u > u 3 , (68) 

corresponds to the scattering solution discussed in subsection C where u = Cf/B 2 approaches 
infinity well before or after the collision, while it reaches its minimum value 113 at some instance 
during collision. In contrast, the case 

27 

7 > 7c = y, Ui < u < u 2 , (69) 

may lead to a motion where u oscillates between the values u\ and u 2 . Also taking into account 
Eq. (|3"Tj) we conclude that all three C\ , C 2 and C 3 would then remain bounded. The aim of this 
subsection is to explicitly construct a special class of solutions that realize Eqs. (|69p and thus 
lead to bounded (quasiperiodic) three-vortex motion. In order to achieve the two conditions 
displayed in Eqs. ((69]) we consider an initial configuration where the three vortices (1,2,3) lie on 
a straight line in the order indicated and thus the relative distances are such that C\ + C 3 = C2 . 
Specifically, we choose 

d = (2 + e)b, C 2 = (3 + e)6, C 3 = b (70) 
and initially place the three vortices along the negative x-axis at points such that 

xi = —C% — C 3 , x 2 — -C 2 , x 3 = —C 3 , 

2/1=2/2=2/3 = 0. (71) 



Then the origin of the coordinate system coincides with the guiding center defined by Eqs. 
Our aim is to calculate the vortex trajectories emanating from the initial configuration defined 
by Eqs. (jTU]) and ([7T |) . 

We first consider the two conserved quantities calculated from such an initial configuration, 
B 2 = ^=(fy£) 2 & 2 , T = -Cf + Ci + C! = 2(3 + e)b 2 , (72) 

and thus 

r 2(3 + £ ) 3 

1= B~^l2Tef (73) 

is now the important dimensionless quantity that controls the behavior of the roots of the cu- 
bic equation. Note that e = leads to the critical parameter j c — 27/2 of Eq. (|67|) while, 
interestingly, 7 > 7 C for both positive and negative values of the parameter e. 

A numerical solution of the initial value problem reveals an interesting picture. For e > 0, we 
encouter a situation of the type described by Eq. (|6"8"|) which corresponds to a scattering solution 
where the vortex-antivortex pair (1,2) moves off to infinity leaving behind vortex 3 which relocates 
to a new final position. However, e < leads to a situation of the type anticipated by Eq. (|6"9")) . 
The resulting bounded (quasiperiodic) three- vortex motion is illustrated in Fig. [S] for e = —1/4. 
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FIG. 6: Trajectories of three vortices in bounded (quasiperiodic) motion emanating from an initial 
configuration specified by Eqs. (|70[) and (|7ip with e= — j. 



Although the present calculation is performed in an infinite magnetic film, this type of bounded 
motion may be relevant for the dynamics in finite, disc-shaped ferromagnetic elements. 

In the special limit e — > the three vortices lie on a straight line for all times at constant 
relative distances: 

d = 2b, C 2 = 36, C 3 = b. (74) 

Such a linear configuration rotates rigidly around a fixed guiding center with constant angular 
frequency 

_ Mi _ d9 2 _ d9 3 1 
W ~ It ~ dt ~ dt ~ 6& 2 ' 

a relation obtained by applying Eqs. I|36[) to the present degenerate case. 



V. CONCLUSIONS 



We have studied the equations of motion for magnetic vortices which are a ppr oximated as 
point vortices. These are a generalization of the Helmholtz-Kirchhoff equations [id . fill derived 
for ordinary fluid vortices. They differ from them in that magnetic vortices have an additional 
characteristic, their polarity, which enters in the equations. The main body of the paper is 
devoted to the case of three magnetic vortices which constitute an integrable system. For the 
integration of this system we follow the steps of the seminal paper of Grobli [12j modified to 
include the vortex polarity. 

Magnetic vortices are commonly observed in magnetic elements with typical dimensions of the 
order of a few micrometers or hundreds of nanometers. A single vortex spontaneously created 
at the center of an element can be driven by a magnetic field or an electrical current to produce 
excitations in its vicinity. It has been observed in various numerical simulations that such 
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excitations often have the form of a lump of magnetization opposite to that of the original 
vortex. In some cases the lumps are amplified to become vortex-antivortex pairs with clearly 
distinct constituent vortices. Such a process was apparently observed in experiments where 
the original vortex was driven by a magnetic field [1,13] or an electrical current [j|. The fact 
that the created vortex-antivortex pair has opposite polarity to the original vortex should most 
probably be attributed to the effect of the magnetostatic (dipole-dipole) interaction which is the 
only magnetic interaction that primarily favors domains of opposite magnetization. From these 
examples it follows that the three- vortex configuration, given in Eq. (|27[) . on which we have 
mainly focused in this paper, is a particularly interesting one for the purposes of understanding 
vortex dynamics in magnetic elements. 

In view of the above comments it is interesting to compare our results with simulations of 
the Landau-Lifshitz equation. We have used the Landau-Lifshitz equation with an exchange 
interaction and an on-site easy-plane anisotropy so that vortices are the relevant excitations of 
the system. Some results have been published in Refs. [1,0]. For the case of vortex-antivortex 
pairs where the vortices are well-separated the simulations show very good agreement with the 
analytical solutions for the vortex trajectories. It is worth discussing separately the special 
case of three-vortex collapse of Subsection IVB which would appear to be specific to the point 
vortex system. We have simulated the case of zero angular momentum within the Landau- 
Lifshitz equation and have observed vortex trajectories very similar to those shown in Fig. [3] 
The original vortex and the antivortex are annihilated when the three-vortex system collapses 
to the origin and an isolated vortex 1 is the final outcome of the simulation. The results of the 
simulation are closer to the analytical solutions when the vortex-antivortex pair is large, i.e., the 
vortex and the antivortex are initially well separated. 
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